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Abstract. Global existence of weak solutions to the Navier-Stokes equation 
in a cylindrical domain under the slip boundary conditions and with inflow 
and outflow was proved. To prove the energy estimate, crucial for the proof, 
we use the Hopf function. This makes us possible to derive such estimate that 
the inflow and outflow must not vanish as t — ¥ oo . The proof requires estimates 
in weighted Sobolev spaces for solutions to the Poisson equation. Finally, the 
paper is the first step to prove the existence of global regular special solutions 
to the Navier-Stokes equations with inflow and outflow. 



1. Introduction 

We consider viscous incompressible fluid motion in a finite cylinder with large 
inflow and outflow, assuming boundary slip conditions. Hence, the following initial 
boundary value problem is examined. 

v t +v-Vv-dxvY(v,p) = / in O T = f2x(0,T), 

div v = in Q T , 

v ■ n — on Sf , 

(1.1) vn ■ B(v) ■ f a + yu ■ f a = 0, a = 1,2, on Sf, 



v ■ n = d on S 2 , 



tT 



n-B(v)-f a =Q, a =1,2, on S£ , 

V \t=Q = W ( ) hl fi > 

where O C R 3 is a cylindrical domain, S = dfl, v is the velocity of the fluid 
motion with v(x, t) — t), V2(x, t), v^(x, t)) E M. 3 , p = p(x, t) 6 K. 1 denotes the 

pressure, / = f(x.t) = (fi(x, f), fi(x, t), fo(x, t)) G 1R 3 - the external force field, 
x — (x±, X2,xs) are the Cartesian coordinates, n is the unit outward vector normal 
to the boundary S and f a , a = 1,2, are tangent vectors to S and • denotes the 
scalar product in M 3 . We define the stress tensor T(v,p) as 

T(v,p) = vB(v)-pI, 

where v is the constant viscosity coefficient and I is the unit matrix. Next, 7 > 
is the slip coefficient and B(t>) denotes the dilatation tensor of the form 

B(v) = {Vi, x . +« J > j }i,j=l,2,3- 
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We assume that O C i 3 is a cylindrical type domain parallel to the axis X3 with 
arbitrary cross section. We set S = Si U S 2 where Si is the part of the boundary 
which is parallel to the axis X3 and S 2 is perpendicular to X3. Hence 

Si = {x G M 3 : <po{xi,x 2 ) = c , ~a < x 3 < a}, 
S 2 (-a) = {x G M 3 : (p (xi,x 2 ) < c , x 3 = -a}, 
S 2 {a) = {x G K 3 : <Po{xi,x 2 ) < c , x 3 = a} 

where a, Co are positive given numbers and fo(xi,x 2 ) = Co describes a sufhciently 
smooth closed curve in the plane X3 — const. 




—\ — 




Figure 1. Domain Q,. 



To describe inflow and outflow we define 

di = -v- n\s 2 (-a) 
d 2 =v- n\ S2{a) 

with di > 0, i = 1, 2. We infer compatibility conditions 



(1.2) 



(1.3) 



didS 2 = / d 2 dS 2 - 

S 2 (-a) Js 2 (a) 



The aim of this paper is to prove the existence of global weak solutions to problem 
(jl.ip without restrictions on magnitudes of external force /, initial data v(0), inflow 
di and outflow d 2 . We would like to show the existence of such solutions that the 
flux does not have to vanish as t — > 00. The presented in our paper method would 
allow us to prove the existence of global regular solutions in the cylinder (in the 
meaning of |RZ3j ) which are much more general than in |Klj . |K2] . [Z] because in 
these papers the flux must converge to zero sufficiently fast. 

We define a space natural for the study of the weak solutions to the Navier-Stokes 
equations: 



v 2 °(n T ) = {u : 



lv 2 °(fi T ) 



eSS SUp ||u||L 2 (f2) 

te(o,T) 




1/2 



^ u \\l 2 (n) dt 



< 00}. 
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To simplify the notation, we do not distinguish between norms of scalar and 
vector function and we write 

3 

ll/H := Y, ll/ill fOTan y /=(/i,/a,/s). 

We also use 

:= 11*11 + 11*11 

for inflow di and outflow d 2 . 

Theorem 1. Assume the compatibility condition 11. 3\) . Assume thatv(0) G L 2 (fl), 
/ € L 2 (0,T;L 6/5 (O)), d, G £oo(0, T; Wp~ 1 ^ p (S 2 )) f] L 2 (0, T; W 2 1 ^ 2 (S , 2 )), § + ± < 

s,p > 3 or p = 3, s > |, di,t G £2(0,T;Wg^g(jS2)),i = 1>2< T/ien t/iere exists a 
weak solution v to problem 1 1.1]) such that v is weakly continuous with respect to t 
in L 2 (fl) norm and v converges to vo as t — >■ strongly in L (ft) norm. Moreover, 
v E V2°(fi T ), v ■ ? a £ £2(0, T; La(Si)), a = 1, 2, and u satisfies 

2 

IMIvJ(n«) +7$I / Tallin < 2||/||| 2(0;t . i6/5( n)) 
(1.4) «=l ° 

+ <p (sup||d|| wrVP(S2) ) (\\d\\\ 2[0 t . w ^ {S2)) + IKH 2 i2(0>t;</ /B (g2)) ) + lk(0)|li 2( o) 

where tp is a nonlinear positive increasing function of its argument and t < T. 

Theorem 2. Assume the compatibility condition \1.3]) . Let f 6 L 2 (kT, (k + 1)T; L 6 / 5 (Q)), 
di e ioo(i?+;W^" 1/p (5 2 )) ni 2 (fcT,(fc + l)T;V^ 2 1/2 (5 2 )), where | + § < s,p > 3 
or p = 3, s > |, and d^ € L 2 (kT, (k + 1)T; Wg/g (5a)), £ = 1, 2. Let its assume that 

\H0)\\L 2 (n)<A 

for some constant A and 

r (k+l)T / x Ak+1)T / s 

; 7 . n/iii 6/5 (n) + * ^ P iidii wrl /, (S3) J y feT ^n^/ a{&) + ii rf *ii<^( S2)/ 

< (1 - e~ vT )A 2 

for k € No, where tp is a nonlinear positive increasing function of its argument. 
Then there exists a global weak solution v to 1 1.1]) such that 

v e V 2 °(fl x (kT, (k + 1)T)) Vfc € N = NU {0}, 

and 

IMI* ?( nx<kT l t)) < 2 / ll/HL /5( n)^ + ^ 2 

(L5) ( \ r ( \ 

+^ ^uplldH^-v,^) ) (JMI^/^ + \\dt\\ w y /!iS2) ) dr 
forte (kT, [k + l)T]. 

The main step in this paper is estimate (|2.15[) - see Lemma 2.2. To derive it, 
we use the Hopf function (see [E], [G]) and estimates in weighted Sobolev spaces 
(see [RZ1J, RZ2].) The estimate is such that we can show global estimate (|4.4[) and 
prove global existence without assumption of vanishing of the inflow-outflow and 
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the external force. The paper makes possible to generalize the result from [RZ3 
into the inflow-outflow case. 



2. Estimates 

To show the existence theorem, we need to obtain the energy type estimate 
and for this purpose, we have to make the Neumann boundary condition (| 1.1 [ Ik 
homogeneous. 

To this end, we extend functions corresponding to inflow and outflow so that 
(2.1) d l \s 2 (a i ) = di, i = 1,2, ai = -a, a 2 = a 

We introduce the function ij, see [L]. 

< a < pe- 1 ' 6 ee r, 




r]{a;e,p)={ -eln~ r < o < p, 

p < a < oo. 



We calculate 

!0 < a < r, 

-- r<a<p, 
p < a < oo. 

£ 

so that \r) (a;e,p)\ < —. We define functions rji on the neighborhood of S2 (inside 
a 

n)-. 

rji = ■q{a l ]£,p), i = 1,2, 
where ai denote local coordinates defined on small neighborhood of S^cti) : 

a 1 = a + x 3 , a- 2 = a~ x 3 

and we set 

2 

(2.2) « = E d ^' 

b = ae3, e 3 = (0,0, 1). 
We construct function u so that 

(2.3) u = v - b. 
Therefore, 

divu = — div & = — a X3 in f2, 
u ■ h — on S. 

Then, the boundary condition for u is homogeneous. The compatibility condition 
takes the form 

/ a, X3 dx = - a\ X3 =- a dS 2 + / a\ X3=a dS 2 = 

Jn Js 2 {-a) Js 2 {a) 
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We define function ip as a solution to the Neumann problem 

A.tp = — div b in £1, 
n ■ V95 = on 5, 



(2.4) 



tpdx = 0. 



Next, we set 

(2.5) w = u — V<p = v — (b + Vyj) = v — S. 

Consequently, (w,p) is a solution to the following problem 
Wt + w ■ Vw + w ■ + S ■ Vw — divT(w,p) 

= f -S t -S-VS+ udiv B(S) = F(S, t) in f2 T , 

div w = in f2 T , 
w ■ n = on S ,T , 
vn ■ 0(w) ■ T a + jw ■ T a 
= —vh ■ B(<5) ■ f a — 70" • ? a = Bi a (S), a = 1, 2, on Sf, 
n-B(w)-f a = -n-B(6)-f a = B 2 a(S), a =1,2, on 5 2 T , 
w\ t=Q = v(0) - 6(0) = w(0) in ft, 

where div S — 0. Moreover, we set 

n\ Sl = n\ Sl = { -^.^f\ f 2 \ Sl = (0,0,1) = e- 3 , 

s 2 (-a) = ~e 3 , n|s 2 ( Q ) = e 3 , n|s 2 = e~i, f 2 \s 2 = e 2 



where ei = (l,0,0),e 2 = (0,1,0). 
We define a weak solution to the problem 

Definition 2.1. We call w a weak solution to problem \2. 6}) if for any sufficiently 
smooth function ip such that 

div^|o=0, ip-n\ s = 

the integral equality 

/ Wt ■ ipdxdt + / H(w) ■ ipdxdt + v I H)(v) ■ H)(ip)dxdt + 7 / w ■ f a ip ■ ? a dS\dt 

2 

— / / B aa i\) ■ ? a dS a dt = F ■ %\)dxdt 

holds, where 

H[w) = w ■ \7w + w ■ V(5 + S ■ Vw. 

Lemma 2.2. Assume the compatibility condition \1.S\) . Assume that f G L 2 (0, T\ L 6 / 5 (fl)), 
di € Loo(0,r;Wp S ~ 1/p (S 2 )) n L 2 (0,T;^ 2 1/2 (S 2 )), where | + § < s,p > 3 or 
p = 3,s > §, dj |t € ^(O,^^ 1 ^ 6 ^)),? = 1,2, w(0) € L 2 (fi). T/ien /or a weafc 
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solution to &2.6]) holds 
2 ,t 



"•||?.:v.iM+7S / H u, - T «lli 3 (s 1 ) < 2 ll/ll! 2 (o,t ; L 6/5 (fi)) 

a=l 170 



2 

n 

(2.7) 

where t < T,d = (dijtfc) and ip is a nonlinear positive increasing function of its 
argument. 

Proof. We use ?/> = w as a test function in a weak solution definition and apply 
definition of F to obtain 

+ f (w ■ V6 ■ w + 5 ■ Vw ■ w)dx — [ div T(w + S,p) ■ wdx 



2 "MM") 



if — St — S' V<5) • wdx 



We use boundary conditions (|f .11) on Si and on 52 to reformulate the third 
integral on the l.h.s. of the above inequality as follows 

div T(w + S,p) ■ wdx = / div [^IB)(w + S) — pi] ■ wdx 
Jn 

= / div [i/lD)(w + §)] ■ wdx — / p • Vw = / Dij(w + S)wj :Xi dx 
Jn Jn Jn 

Dij(w)w jtXi dx + f Dij(5)wj >Xi dx = - f \Dij(w)\ 2 dx+ f D ij {5)w j _ Xi dx 
n Jn £ Jn Jn 

Then, we apply Korn inequality to have the estimate 

Id 2 

9 ;dMli a (n) + "IMl!ri(n) + 7 H w ' f «lli 2 (s 1 ) 



2 



(2%_ /" ( w .V5-w + 6-V w -w)dx + cY\\5-f a \\l 2iSi} +cmS)\\ 2 L 2i n) 

Jn a=l 

+ [ (f-5 t -5-VS)wdx 
Jn 

Next, we focus on the integral 

/ 5 ■ Vw • wdx = I (b + Vy>) ■ Vw • wdx — b ■ Vw ■ wdx + / V<p • Vw • wdx 
Jo Jn Jn Jn 

= h+h 

We can estimate I\ by Holder inequality and definition of b 

\h\ < ||Vw|| £2(fi) ||u;|| L6( n)||6|U 3 (n) < 4w\\m(n) ll fe lli 3 (s 2 ( P )) 
< cp 1/6 \n\ 2 H H n)\\b\\ L6{ s 2{p )) < cp 1/6 |H| 2 ffl( n)ll<5|U 6 (n) < cp 1 /«|j lt ,||| rl(n) ||d||^ 1(n) 
where 

S 2 (p) = {x G : x 3 £ (-a, -a + p) U (a - p, a)} = S 2 (p, fli) U 5 2 (p, a 2 ). 
We estimate i 2 as follows 

(2-9) J a = / V^.Vto-«;da!<||Vfp|U,(n)||ti;|Ue(n)l|V«;|| £a (n) 
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where 



l|V</?|U 3 (o) < c\\Vtp\\ L3 _^ m < cllV^V^II^^^n) < c\W\\lI :1 _^,(Vl) 

< c||div&||i (tl ) 



and we denote 



\m\Li 



„(n) = (J2 I \D:u\ 1 \jim =h2 \(dist(x : S 2 (a l ))rdx) 1 ^^eR,pe (l.oo). 

\a\=k J 

To estimate the last norm, we have used the result of (RZ2 on Poisson equation in 

2 

weighted Sobolev spaces and choose - < 1 — // < 1. With /i = 1 — // we have 

E /. n 3 ^h* ) + (E /~ i4x 3 i 3 ip(^)i 3 ^x 



<cVe(sup/ \di\ 3 dx' I \dcrA + V ( sup / \d itX3 \ 3 dx'[ a^dxn 

i=l \ x * JSi{ ai ) Jr °t ) 7^1 \ X3 • , S 2 {a i ) ' 



1/3 



1/3 



/n 



< cep" 2 / 3 sup ||d|| i3(S2) + cp^ +1 / 3 sup ||d, X3 || L3( s 2 ) 

213 ^3 



where oi — dist{S2(ai),x},x £ S2{a%,p)- We note, that the last bound holds for 
jj, > | since for p = | the r.h.s. takes the form 



csup||J|| i3(S2) +cp sup \\d, X3 || L3 (s 2 ), 

a:3 ^3 

which can not be made small for large d. Then, 

ep^- 2 / 3 sup || J|| i3 (S 2 ) + P M+V3 sup HJ,^ || ia(Ss) 



h < c 



IHIffi(n) 



Next, we consider the term 



(w ■ V<5 • w)dx — (w ■ Vb ■ w)dx + J (w • Wip • w)dx = I3 + 14. 
n Jn J q 



For I4, we have 



div (w • V</3 • w)(ia; — (w ■ Vw ■ Vf )dx 
n Jn 



\h\< 



< / \n- V(p -w z \dS + / \Vtp ■ (w ■ Vw)\dx < / \V (p ■ (w V w)\dx 
Js Jn Jn 

so I4 can be treated in the same way as I2 and therefore 



(2.10) |J 4 | < c 



ep^ 2 ' 3 sup ||<f|| L , (a) + P M+1/3 sup \\d 



)X 3 Ni 3 (S 2 ) 



Lff!(o)- 
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On the other hand, using b — ae~3 = ^Zj=i <^i e ~3j we find the bound for Is 



\h\ < \y~] I w ■^{d % r) i )w z dx\ 

J S 2 (p,a;) 



< IE 
* El 

i=l 
2 



(u> • X7dir/iW3 + w ■ VriidiW^)dx\ 



S 2 (p,ai) 



|w • Vdi?7i||it;3|c?x - 





W 3 ~ 
— W 3 di 


d(Jidx\dxi 


'S 2 (p,ai) 







S 2 (p,a,) 

i6(S 2 (p,a,)) ll ' i ' ;i|l L3(S2(p,a l )) 11 v "» 11 ia(Sa(p,o 4 )) 



< cJ^IHI^^^^Hiflsll,..,?.,.. J|Vd 

i=l 



i6(S , 2 (p,a i ))ll u ' J llL3(S2(p,a i )) 



2 



dx\dx2 / dcr. 



S 2 (p,a;) 



7«3 



< cp^^iHi? II V(/. 



i=l 
2 



i 6 (S 2 (p,a,)) H ,|| i2(S 2 (p,a l )) 

+CeEH U 'lli 6 (S 2 (p,a0)ll Vu '3llL 2 (S 2 (p,a0)ll^llL3(S 2 (p,a 1 )) 
i=l 

< C(p 1/6 +£)||w||^ 1(n) ||J|| w i (f2) . 



Thus, we can summarize estimates for /i — I4 to conclude that nonlinear term 
(HH) is bounded by 



(w • V<5 • u> + (5 • Vu> • w)dx 

2 

(2-11) < c|HI^) f^- 2 / 3 sup ||rf|| i3 (s 2 ) + P M+V3 sup || J,. 3 \\ L3{S2 

\ X3 X3 

+ (p 1/6 +e)\\d\\ P ^\\d\\ 
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Next, we examine the second term on the r.h.s. of (12.8 



Eli*- 

a=l 



< Edl^lli^o + IIVvr«||i a(Sl) ) 



< 



ll a lli a (sx) +c ll v *'llwj 

2 



3/2^) 



^ E 11*111,(50 + c ll div& lli 3/2 (o) 

1=1 

2 

^ c ll<*llw* a (n) + C E H V (*^)llL/ 2 (n) 

i— 1 
2 

^ c ll^llw 3 i /2 (o) + C E (ll v *^lli 3 /2(n) + ll*V»»||i a/a( n) 



i=l 
2 



< c||d||^i (n) +c^||J i V??i|| 

2—1 



i 3 /2(n)- 



The last expression we calculate in details: 



En* Vr * 



i|lL 3/2 (f2) 



f—a+P r 


dx 


3/2\ 


1 dxz I dx 






l-a+r JS 2 (a 1 ) 


a + x 3 





4/3 





d 2 


3/2\ 


/ dx$ / dx' 






wa-p J 82(0,2) 


a - X3 





4/3' 



su pIMiIIL/2(s 2 («o) 



1 



a + x 3 



< C£ 2 snp||d||i 3/2(S2 



+ SU P ||J 2 ||| 3 /2 (S 2 (a 2 )) 

2:3 Wa-p 

dy \ 4 ^ 3 2 .. 7.. 2 

^372-J ™ Pl|d|li 3/ 2 (S 2 ) 



1 


3/2 \ 


dX3 J 


a - x 3 





3/2 \ 4 /3 

dx 3 

4/3' 



4/3 



1 1 



-1/2 l/2 



< C£ 2 sup||d||| (S2 _ 

X-3 P ' 



1 [e l/ 2e _ 1]4 /3 < c _|L e 2/3 £sup ||J| 



^n-uL 3/2 (S 2 )- 



Combining inequalities above, we infer 



2 



^\\S- f a \\i 2{Si) < c\\d\\U m + c-^e^ sup \\d\H {S2) 

a=l P X3 
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We estimate also the term 

wmwUn) < \\m\\l 2m + mv<p)\\i 2m 



< 



< 



L 2( a) + IM>^ll! 2( n) ) + l|VVlli 2( n) 



i=l 
2 

E (ll v *»»lli a (n) + ll*VTji||i a(n) ) + ||div6||i a(n) 



< £^2[\\Vdir}i\\ L2in) 

i=l 
2 

< c2||*||^ 1(n) +e 2 c 



/ dx 3 / dx' 

' a p J 82(0.2) 



\diVrn 



\L 2 (tt) 



r-a+P r 




/ dx$ / dx' 




1 —a+r J Si (a\) 


a + x 3 


2 





a - x 3 



; E 

1=1 

2 

; E 



< c 

1=1 

2 

c 
1=1 

2 

* c E 

1=1 

2 

* c E 

i=i 



illWf(O) su Pll^lll 2 (S 2 ) 



ay 

,2 



Kllvt^n) +£ 2 nup||di||i s 



(s 2 ) 



p 

V 

1 1 



l*llwj(n) +^ s up||d l |li 2( 5 2 )^(e 1/£ - 1) 



e 2 

\ d i\\wu Q) + -e 1/£ sup ||rf*||i a(Sa) 



Analyzing the last integral on the r.h.s. of (|2.8p we have 

jf (/ W)uKfc < ei||t«||i 6(n) + c(l/ £ i)(||/||! 6/5(n) + ||^||| 6/5(n) ) 



S ■ V<5 • wda; 



We estimate ||5t||i 6/5 (n) as follows 

INU 6/5 (n) = \\bt + Vip t \\ L6/5m < ||dt|U 6/5 (o) + ||div6 t || L6/5(n) 
< ||*IU a/ ,(n) + l|Vd t || L6/5( n) + ||*V»7||£ a/5 (n) 

( rP dx 3 

6 

It Xq 



5/6 



< l|dt||wj /8 (n) +£Sup||rft|U 6/5 (s 2 ) ( / 

< l|dt|l^ /5 (^) +£-r^e 1/6e sup||J t || i6/5(S2) 
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since 

5 / 6 . .. .. s 5/6 



It X 3 

Finally, we examine 



J r xV 5 W /8 P 1/5 J P 1/6 ^ ) 



< 



l|V5|U 2 (n)lhlli 6 (n)ll*IU 3 (n) < MHIi 8 (n) + c (V £ 2)ll<*llw J 1 (fi) 
<ealltull? f lldlli,,,^ + ie 2 / £ suDllJll1 



£ 2 |klli 6( n) + c(l/e 2 ) l|d||^(o) + 72 e /£ SU P Nll 2 (s 2 ) 

\ P x 3 

We summarize above estimates to rewrite (|2.8j) as follows 
Id 



2^IIHII 2 (n) + f|Mliri ( n) + 7^3 Ik • ^'Wo 



< \\ w \\m(n) 



2 



2/3 sup||d|| L3( s 2) + // 1+1/3 sup||d, :C3 || ia(Sa) 

X 3 X 3 

+ (p 1/6 +e)\\d\\ ff!(n) + £i + £2 

(2 - 12) + ll/lli 6/5( a) + NlL ( n) + l|d||^(n) + lld|lwj ( n) 

+ HVd1i 6/6(n) + ||Vd|| 4 L2(0) + ||d||^ /a(n) + 11*11^(0) 

2 4 

+ e ~z lte sup ||d~||| 2(S2) + ^e 2 / £ snp ||d~||i (S2) 

2 2 

+ 47J e2/3£su Pll J HL /2 (s 2 ) + TT7I el/3esu Pll J *l1l/5(S 2 ) 

We apply Sobolev anisotropic imbedding (see I BIN) , Ch.3, Section 10) to estimate 
sup X3 ||d||i 3 (5 2 ) and sup X3 ||d )X3 ||l 3 (s 2 ) with some W£ norm and calculate 



2 77 - + + - < 1 for P>3 

\p 6 J s p s s 

2( 1 -1) 1 - + 1 - 1 - + 1 -<1 for P = 3 



Then, 

3 1 4 

(2.13) - + -<s for p > 3 or p = 3, s > -. 

p 3 " 3 

We set /i > |, then since p < 1, we observe that p^ +1 / 3 < p 1 / 6 . Then 

£ ^- 2 / 3 sup||d|| i3(S2) +^+ 1 / 3 sup||W|| i3(S2) + (p 1/6 + E)||dl|w 1(n) +p 1/6 ||d|| iflcn) 

< {ep^ + P ^ + 2p^ + E)\\d\\ w;(n) < (2e + 3p^)\\d\\ W s (n) 
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We put 

V 

£ = = . 

!5Hw2(n) 

1/6 _ v 



(2.14) p 1/e = 



l5\\d\\ W s {n) 



£l+£-2 = -x, 
6 



with p, s satisfying (|2 . 13|) . Therefore, 

ep^ 2 ' 3 sup \\d\\ L3(S2) + p^ 3 sup ||Vd|| £s(Sa) 

X3 X 3 

+ (P 1/6 + e)\\d\\ w i m + P 1/e \\d\\m(n) + e l +e 2 <^ 
and formula (|2.12[) assumes the form 



d 



Q — 1 

< 2 ll/lll 6/5 (n) + V(\\d\\wun))(\\d\\ 2 wjin) + H*H^ /B (n)) 
+<p{\\d\\wi[a)) ( sup \\d\\l 2{S2 ) + SU P Klli 6/B (s 2 ) 



where ip is a nonlinear positive increasing function of its argument. We use Sobolev 
imbedding 

sup||J|| L2( 5 2) < c||d||wi( n ), 
suplKIU (S 2 ) < c\\d t \\ w i m 

and hence 

d 2 

-dM!z, 2 (fi) +"\\ w \\m(n) +7^2\\w-f a \\l 2(Sl ) 

(2.15) at ' a =i 

< 2||/||i 6/5( o) + ¥>(||d|k;(n)) (\\d\\w im + IWI^ /5 (n)) 
Integrating (|2.15|) with respect to time we obtain 

NlSj(n«)+7$^ / W W ' f «lli a (Si) d * ^ 2 H/Hi 2 (0,t;L 6/B (fi)) 

(2.16) «=i Jo 

+ <p(sup\\d\\ W s {n) ) (llrflll^o.t^Kn)) + \\dt\\ 2 L2{Qtt , w i /&(n)) ) + 11^(0)111,(0), 

where ^ + ^<s,p>3orp — 3, s > |. 



□ 



In this section, we use the Galerkin method to prove the existence of weak 
solutions to the problem (|2.6[) . We follow ideas from [L], chapter 6, section 7. 



3. Weak solutions to (12.6 
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Namely, we introduce the sequence of approximating functions wn given as 

JV 



(x,t) = Y,C kN (t)a k (x), 



fc=i 



where {a k }^Li is the system of orthogonal functions in L 2 {VL) f] J" (12). Here, J® (12) = 
{/ € if 1 (12) : div/ = 0} and {a fc }^L 1 is the fundamental system in if 1 (12) with 
sup^gf; |a fe (x)| < oOjSup^g^Q ^(^Ol < °°- The coefficients CfcAr(O) are defined by 

C fc jv|t=o = {w ,ak), k = l,...,N, 
and the function w N satisfy the following system with test functions a k : 
'1 d 



w N a k + w N ■ Vw N a k + 6 • Vw' v • a k + w ly -\7S-a k + M{w ly )B(a k ) 



..N k 



.N k 



N 



,.N\ 



2 dt 



> f w N ■f j a k f j dSi 1 = ( V / B aj a k ■ fj dS a + f 
JSi J \ jt<T=1 Js a Jn 



F ■ a k dx 



for k — 1, . . . , N. Then, w would be the weak solution to (|2.6[) . 

With (/, g) = J n fgdx and {f,g)s — Jg fgdS this can be rewritten as: 

(w? , a k ) + (/ • Vw N , a k ) + {5 ■ Vw N ,a k ) + (w N ■ VS, a k ) 



+v{0{w N ), 0{a k )) + 1 (w N ■ fj, a k ■ fj) Sl 

2 



]T {B a] ,a k -T 3 )s„+{F,a k 

c,J=l 



, k = l,...,N. 



Thus, 



(3.1) 



dt 



w N , a k )+ (w N • W v , a k ) + (5 ■ Vu/ V , a k ) + {w N ■ Vtf, a k ) 



+ ),B(a fc )) + iK ' r 3 ,a k ■ f,) Sl 



= ( B °i> ak ■ f l)s„ + (F, a k ), k = 1,...,N. 

The above equations are in fact a system of ordinary differential equations for the 
functions CkN{t)- The properties of the sequence a k imply 

N 



(^t)la,n = E C *ArC*)' 



k=l 



On the other hand, we can obtain the a priori bounds for the approximative solu- 
tions w N of the same form as (|2.16p : 

\w N \ 2 V 0W)= SU P \w N kn+ [ \Vw N \ 2 ,ndt 

2 y ' 0<t<T Jo 



(3.2) 



< 



H/lliWi» + <P( SU P Mlw'fli)) / (l|d|lwri(n) + 11*11^,(0) 



- 8 /s(n) 

2 



0<t<T 



dt 



\W N (0)\\ 2 T .^ < C, 
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where | + | < s,p > 3 or p — 3,s > |. Therefore, sup 0<t<T |CfcAr(£)| is bounded 
on [0, T] and w N are well defined for all times £. 

Let us define now if> N ,k = (w N (x,t),a k (x)). This sequence is uniformly bounded 
by (|3.2[) . We can also show that it is equicontinuous. Namely, we integrate (|3.ip 
with respect to £ from £ to £ + A£ to obtain 

ft+At 

|Vjv, fe (£ + At) - ^jv,fc(t)| < sup |a fc (a;)| / (\w N ■ Vw N \ 2 ,n + \6 ■ Vw N \ 2 , n 



ft+At 

\w N ■ V<y| 2 ,n + \F\ 2 ,a) dt + v\Va k \ 2 , n J \Vw N \ 2t ndt 

f t+At / 2 \ 

7sup|a fe (x)| J \w N •f j | 2 , Sl + ^ |S CTJ | 2 ,s„ ] d< 



| to 

J>=1 

,.N\ /iv7„„Ai i iv7Jti ^ i „„„ ixi ir7„..JVi 



dt 



< sup \a k (x)\VAt sup \w N | 2 , n (|Vto w | 2>fr + |V<5| 2 , QT ) + sup \5\ 2M \W \ 2 ^ T 

ft+At 

+ sup|a fe (a;)| / \F\ 2jQ dt + p\Va k \ 2 , n VM\Vw N \ 2i nT 

i£!I Jt 

( ft+At 2 \ 

+ 7 sup|a fe (x)| I y/At\Vw N \ 2 , nT + jf X>;k S )j 

(/•t+At 2 

We can see that for given k and N > k the r.h.s. tends to zero as A£ — > uniformly 
in TV. Thus, it is possible to choose a subsequence N m such that ^N m ,k converges 
with m — > oo uniformly to some continuous function ^ for any given k. Since the 
limit function to is defined as 

oo 

w(x,t) =^2i/j k {t)a k (x), 

k=l 

then we conclude that (w Nm — w,ip(x)) tends to zero as m — > oo uniformly with 
respect to £ G [0, T] for any ip £ J 2 (&) and w(x,t) is continuous in £ in weak 
topology. Moreover, estimates (|3 . 2[) remain true for the limit function w. 

We will show that {w Nm } converges strongly in L 2 (fi T ). To this end, we need 
to apply the following version of the Friedrichs lemma: for any e > 0, there exists 
such N £ that for any u £ W 2 (0) the following inequality holds: 



2 

l2.fi- 

k=l 



\u\\l n <Y,(^a k )+s\\Vu\ 



This in terms of u — w Nm — w Nl reads 

\\w N ™-w N % nT < V / (w N ™-w Nl ,a k )dt + e\\Vw N ™-Vw N i\\l nT . 
By (O, we have 

\\Ww Nm -Ww N '\\ 2 2 nT < 2C 2 



LARGE FLUX FOR NSE 



15 



for some constant C. The first integral on the r.h.s. for given number N e can 
be arbitrarily small if only m and I are sufficiently large, so it tends to zero as 
to, I — > oo. Therefore, {w Nm } converges strongly in L 2 (Jl T ). 

We summarize the above convergence properties of the sequence {w Nm } : 

(i) {w Nm } — > w strongly in L 2 (n T ) for some w, 

(ii) {w Nm } —> w weakly in L 2 (il) uniformly with respect to t g [0,T], 
(hi) X7{w Nm } Vffl weakly in L 2 {Vt T ). 

With given $ fc = Y^,j=i dj(t)a^(x), the sequence {/-} satisfy the identities: 

/ (i wNm ® k + ( wNm ' ^ wNm + 6 ' ^ wNm + wNm ' v<5 ^ $fc + M ( wNm M® k )^J dx 

2 

+7 / /" • ?j<S> k ■ fjdSo = V / B aj <S> k ■ fjdS a + I F$ k dx. 
J Si a , j=1 Js tr Jn 

Then, we can pass to the limit with m — > oo to obtain the identity for w. 
Conditions divw^ = 0, w N ■ n\sT =0 stay true for the limit function w as well. 

It remains to consider the limit lim t ^o w(x, t). We note, that w Nm satisfy the 
relation (|2.8p (if we use the test function w m ) . This yields 

\w Nm \ 2<n < |itfd|a,n + / (\F\ 2 .n + \B\ 2 ,s)dt. 
Jo 



In the limit m — ¥ oo we obtain 

rt 



\w\ 2 .n < \w \ 2 ,n + I (|-F| 2 ,n + \B\ 2 ,s)dt 
Jo 



which implies 

lim t _i.o|if|2,n < |tt?o|2,n- 
On the other hand, since w Nm tends to w as m — > oo, we have \w Nm — tuob.n - > 0. 
Therefore, \w Nm — wo \ — > weakly in L 2 (f2) as t — > and 

\wo\ 2 ,n < lim t ^ |u>| 2 ,fi. 

We conclude that the limit lim^o |u?|2,n exists and is equal to |u>o I2.fi where the 
convergence is strong - in the norm L 2 (fl). 

Consequently, we have proved the following result. 

Lemma 3.3. Let the assumptions of Lemma \2.2\ be satisfied. Then there exists a 
weak solution w to problem \2. 6\) such that w is weakly continuous with respect to 
t in L 2 (fl) norm and w converges to Wq as t — > strongly in L 2 (fl) norm. 

Since v — w — 5 we conclude the analogous existence result for v formulated in 
Theorem 1. 



4. Global solutions to 
To obtain a global estimate we write (|2.15[) in the form 

^IIHlL(fi) + HHlL(fi) ^ 2 ll/llL/5(^) + ^(Nli^(fi)) (ll<*Uw*(n) + IKHw^fi)) 

where | + i<s,p>3orp = 3,s> |. Hence 



| (|MI W') < 2||/||! 6/5(n) e^ + <p(\\d\\ W s m ) (\\d\\lr m + IKH^ /5 (fi)) e " 
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Integrating with respect to time from ti to ti yields 

IM^IlW' 2 < 2 £ 2 ||/||i 6/5(n)e ^t +11^)11^^^"- 



+V>(BUp||d|| W r ;(n )) ^ (ll^ll^tO) + 11*11*^,(0) 

Thus, 

\\^)\\ m < 2e"^ J** ||/||i 6/B( n)^dt + \\w(h)\\l 2m e~^-^ 

+ V (sup\\d\\ Wi{n) )e-^ £ 2 (||d||^ (n) + Wdtfw^u)) e vt dt 
and this implies 



e vt dt. 



lh(<3)||i a( n)<2/' ta ||/||i. /1 , (n )A+imti)||i a( n)C- , ' ( ' J -'" 

+^(sup||d|| Wp .(n))^' (||d||^i (n) + 11^11^(0) 
Setting ti = and £2 = i € i?+ we obtain the global estimate 

IK*)ll£ 2( n) < 2 / t ||/lli 8/B (n)^+im0)||i 2(n) e-" t 
(4-2) /■* 

Let fc e N. Integrating (|2.15p with respect to time from kT to t E {kT, (k+ 1 )T] we 
get 



dt 



(4.3) 
Therefore, 
(4.4) 

We have also 



My 2 °(fi X (fcT,t)) < 2 / Il/Ili 8/B (n)dr + ||to(*r)Hi a (n) 

J kT 

+<p(swp \\d\\ W s {n) ) J (||d||^i (n) + ll^llw^n) 



Mv 2 °(nx(fcT,t)) < 2 / ll/lli a/6 (n) dT + IH fcr )HL(fi) 
+(^(su P ||d||Hr.(n)) [ {\\d\\wi (n) + l|d t ||^i /6(n) 



■t 

2 ^ r, / 11X112 j_ , ll„./n\ll2 ,-l/T 



e 



2 l 2 (q.) < 2 / ll/Hl ( n)dr + ||w(0)||| 2(n) 
(4.5) Jo t 

+cp(sup\\d\\ W s m ) (||d||^i (n) + lldtllw^ti)) 

We set hi = e~ vT . Let as assume that 

ll«(0)|U 3( n) < A 



dr 
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for some constant A and 

2 f ll/lli. /B(n )^+^(sup||d|| W j (n) ) / (\\d\\w m + \\dt\\w c] m ) dr < (l-e-» T )A 2 

J " J 

Thus, 

\\v(T)\\ L2(n) < A 

so we can control the initial condition for the next time step. This can be repeated 
for intervals (kT, (k + 1)T). Then by (|4.4j) we can prove global existence of weak 
solution such that 

v e V?(Q x (kT, (k + 1)T)) Vfc e N = N U {0}, 

so we conclude Theorem 2. 
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